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Preventing Multipartite Disentanglement by Local Modulations 
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An entangled multipartite system coupled to a zero-temperature bath undergoes rapid disentan- 
glement in many realistic scenarios, due to local, symmetry-breaking, differences in the particle-bath 
couplings. We show that locally controlled perturbations, addressing each particle individually, can 
impose a symmetry, and thus allow the existence of decoherence-free multipartite entangled systems 
in zero-temperature environments. 
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Symmetry is a powerful means of protecting entan- 
gled quantum states against decoherence 0], since it al- 
lows the existence of a decoherence-free subspace or a 
decoherence-free subsystem (DFS) |2|,|ig. In multipartite 
systems, this requires that all particles be perturbed by 
the same environment. In keeping with this requirement, 
quantum communication protocols based on entangled 
two-photon states have been studied under collective de- 
polarization conditions, namely, identical random fluctu- 
ations of the polarization for both photons Q . Entangled 
qubits that reside at the same site or at equivalent sites 
of the system, e.g. atoms in optical lattices, have likewise 
been assumed to undergo identical decoherence 

However, locally-decohering entangled states of two or 
more particles, such that each particle travels along a dif- 
ferent channel or is stored at a different site in the system, 
may break this symmetry. A possible consequence of this 
symmetry breaking is the abrupt "death" of the entan- 
glement [5|. Such systems, composed of particles under- 
going individual or "local" decoherence, do not possess a 
natural DFS and thus present more challenging problems 
insofar as decoherence effects are concerned |6|. Specifi- 
cally, can the multipartite symmetry, broken by "local" 
decoherence, be dynamically restored to allow a DFS? Al- 
ternatively, can we dynamically impose other symmetries 
that protect all particles, not just part of them, from lo- 
cal decoherence? Previous treatments Q have suggested 
to dynamically restore a DFS by subjecting the particles 
to sufficiently frequent, fast and strong pulses, assuming 
that they stroboscopically decouple all particles from the 
bath, irrespective of their coupling strengths. 

Here we address these fundamental questions by devel- 
oping a generalized treatment of multipartite entangled 
states (MES) decaying into zero-temperature baths and 
subject to arbitrary external perturbations whose role is 
to provide dynamical protection from decay and decoher- 
ence. Our treatment applies to any difference between 
the couplings of individual particles to the baths. It 
does not assume the perturbations to be stroboscopic, 
i.e. strong or fast enough, but rather to act concurrently 
with the particle-bath interactions. Our main results are 



to show that by applying local (selective) perturbations 
to multilevel particles, i.e. by addressing each level and 
each particle individually, one can create a decoherence- 
free system of many entangled qubits. Alternatively, one 
may reduce the problem of locally decohering MES to 
that of a single decohering particle, whose dynamical con- 
trol has been thoroughly investigated 0, H, • 

In our treatment we introduce the multipartite system 
and arrive at a general dynamical solution for its deco- 
herence. We then analyze realizations of different sym- 
metries of the decoherence matrix comparing the result- 
ing fidelities of the decohering state. Our system com- 
prises N particles, labelled by j — l,...,N, each having 
a ground state \g)j and Mj excited states, {!«),,} with 
energies, hujj_ n . In particular, nj may enumerate the mo- 
tional states of a trapped ion or atom, 6]. The coupling 
to the zero-temperature bath may differ from one par- 
ticle to another and from one excited level to another. 
For their protection from decay and decoherence, we ap- 
ply a perturbation such that each level is modulated by 
a different AC Stark shift HSj tn (t) and/or by different 
off- resonant perturbing field €j M (t). The total Hamil- 
tonian is the sum of the multipartite system, bath and 
off-diagonal system-bath interaction Hamiltonians |Toj : 

JV Mj 

#0) = h J2 £ k» + *»>(*)] \ n )jM + fiJ2^k\k)(k\ + 

j—1 n—1 k 

N Mj 

''Y.Y.Y.'^ 11 '^*" (\ n )(9\)j \vac)(k\ ® I f7tj + H.cS\l) 
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Here / is the identity operator, fJ-k,j,n 1S the off-diagonal 
coupling coefficient of the n th excited level of particle 
j to an excitation of the bath mode \k) and \vac) is the 
vacuum state of the bath. H.c. are Hermitian conjugates. 

Although our strategy applies in general to any num- 
ber of excitations, simple closed- form solutions are ob- 
tainable for a single initial excitation of the system. Ac- 
cordingly, the complete wave function is: 
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We will denote the first and second RHS terms as the 
bath and system wave function, \^S> B (t)} and |^E' s (i)), re- 
spectively. In order to solve the Schrodinger equation, 
one may eliminate the {(^(i)} amplitudes and trans- 
form to the interaction picture, ending up with an exact 
integro-differential equation. Assuming that these am- 
plitudes are slowly varying on the time-scale of the bath 
response and using the matrix representation, this 
equation has the general solution: 

a(t) = T+e- J ^">a(0) 7 &j „ = e^-"*+ 1 ^ dT ^-^ aj n . 

(3) 

where a = {ctj.n}, T+ is the time-ordering operator and 
J(t) = {Jjf ,nn'{t)} is the dynamically-modified decoher- 
ence matrix, determined by the following convolution: 



Jjj',nn'(t) = 27T / dhjGj tji lTm ' (t-o)K t Jj> ,„ n ' (uj) (4) 
J — oo 

Gjj>, nn '(u>) = IT 2 Hk,j,nt4,j',n> S ( u - ^k) (5) 
fc 

Kt,jj',n»(u) = <j>(w - ^j,n)e t ,j', n '{^ ~ Uj'n') (6) 

etj,»M= / rfre J> (r)e^ T (7) 
Jo 

e J> (t)= e - J - n (t)e-/o^.4r) (g) 

Here Gjj' (w) is the coupling spectrum matrix given 
by nature and K t jji nn i(uj) is the dynamical modula- 
tion matrix, which we design at will to suppress the 
decoherence (cf. Q). This general solution holds for 
dynamically-modified relaxation of a singly- excited MES 
into a zero-temperature bath. 

For example, we may control the decoherence by im- 
pulsive phase modulation, i.e. a sequence of pulsed Stark 
shifts caused by fields whose amplitudes satisfy fj 



(9) 



^ = ' ^ek^,,) _ 1} ^ 10) 

Here [...] denote the integer part, Tj n and 9j n are the 
pulse duration and the phase change for level n of par- 
ticle j, respectively. In the limit of weak pulses, of area 
\0j,n\ < 7T, Eq. (TJU yields e tJ> (w) = e tJ> <5(u; - A i)tl ), 
where Aj in = 8j. n /Tj. n is the effective spectral shift 
caused by the pulses. 

It is expedient to rewrite the fidelity of the evolving 
system, F(t) = |(* s (0)|# s (i)) | 2 , as a product of two 
factors: 



F(t) = F p (t)F c (t) 

N Mj 

F p (t) = \A(t)\ 2 =J2Y,\ a ^W 
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where F c (t) is the autocorrelation function, 

|(\I>' s '(0)|\I/' s '(t))| 2 , normalized by the total excitation 

probability |A(t)| 2 . Thus 1 — F p (t) measures population 

loss from any \n)., whereas 1 — F c (t) is a measure 

of correlation preservation: F c = 1 when the initial 

multipartite correlations are completely preserved. As 

shown below population and correlation preservation can 

be independently controlled in the model of Eqs. (|T|) . i|2")l . 

In the absence of dynamical control, F c (t) decays 

much faster than F p (t) and is much more sensitive to 

the asymmetry between local particle-bath couplings. 

Thus, for initial Bell singlet and triplet states, |^(0)) = 

1 / y/2(\g) A \n) B ± \n) a \g) B ), which do not experience 

cross-decoherence but only different local decoherence 

rates, a A{B) {t) = 1 / ^/2e~ J (t). We find F p (t) = 

( e -2Jx(t) + e -aJfl(t))/ 2 ; F c {t) = (1 + C{t))/2 = 1/2 + 
e-AJW/^ + e -2Aj(t)) j where AJ ( t ) = j a (q _ j B ( t ) 

and C(t) is the concurrence [TTj . 

We shall first deal with N identical qubits, and thus 
ignore the n subscript, i.e. set ujj = loq- We also require 
that at any chosen time t = T, the AC Stark shifts in 
Eq. |J3J satisfy, J Q T dTdj(r) = 2irm, where m = 0,±1, .... 
This requirement ensures that modulations only affect 
the decoherence matrix Q, but do not change the rel- 
ative phases of the entangled qubits when their MES is 
probed or manipulated by logic operations at t = T. 

Without any modulations, decoherence in this scenario 
has no inherent symmetry. Our point is that one can 
symmetrize the decoherence by appropriate modulations. 
The key is that different, "local", phase- locked modu- 
lations applied to the individual particles, according to 
Eq. ©, can be chosen to cause controlled interference 
and/or spectral shifts between the particles' couplings to 
the bath. The Kt jj>(u>) matrices (cf.®) can then satisfy 
2N requirements at all times and be tailored to impose 
the advantageous symmetries described below. By con- 
trast, a "global" (identical) modulation, characterized by 
Ktjj'(uj) = \e t (uj)\ 2 , is not guaranteed to satisfy N ^> 1 
symmetrizing requirements at all times (Fig. ^i). The 
most desirable symmetry is that of identical coupled par- 
ticles (ICP), which would emerge if all the modulated 
particles could acquire the same dynamically modified 
decoherence and cross-decoherence yielding the following 
N x N fully symmetrized decoherence matrix 



J^(t) = r(t) Vj,f. 



(14) 



ICP would then give rise to a (N — l)-dimensional 
decoherence-free subspace: the entire single-excitation 
sector less the totally symmetric entangled state. An 
initial state in this DFS @ would have F(t) = 1 for all 
times, meaning that it would neither lose its population 
nor its initial correlations (or entanglement). 

However, it is generally impossible to ensure this sym- 
metry, since it amounts to satisfying N(N — l)/2 condi- 
tions using N modulating fields. Even if we accidently 
succeed with N particles, the success is not scalable to 
N + 1 or more particles. Moreover, the ability to impose 



FIG. 1: Two 3 - level particles in a cavity, coupled to the 
cavity modes (thin lines) and subject to local control fields 
(thick lines), (a-c) Frequency domain overlap of coupling 
spectrum (dotted) and modulation matrix elements (solid), re- 
sulting in modified decoherence matrix elements (shaded) , for: 
(a) global modulation, (b) cross-decoherence elimination (IIP 
symmetry) and (c) IIT symmetry, (d) IIT symmetry scheme. 

the ICP symmetry by local modulation fails completely if 
not all particles are coupled to all other particles through 
the bath, i.e. if some Gjji (co) elements vanish. 

A more limited symmetry that we may ensure for N 
qubits is that of independent identical particles (HP). 
This symmetry is formed when spectral shifts and/or in- 
terferences imposed by N modulations cause the N dif- 
ferent particles to acquire the same single-particle deco- 
herence r(t) and experience no cross-decoherence. To this 
end, we may choose in Eq. © et,j(w) — e t jS(uj — Aj), 
The spectral shifts Aj can be different enough to couple 
each particle to a different spectral range of bath modes 
so that their cross-coupling vanishes: 

Jjf(t) = j duG(cj)6(u} - Aj)5(u - A r ) -> 0. 

(15) 

Here, the vanishing of Gjf (oj) for some j, j' is not a lim- 
itation. The N single-particle decoherence rates can be 
equated by an appropriate choice of N parameters {A^}: 

jW(t)^\^\ 2 G n (A 3 )^5 n ,r(t), (16) 

where Sjf is Kronecker's delta (Fig. QJ)). The IIP sym- 
metry results in complete correlation preservation, i.e. 
F c {t) = 1, but still permits population loss, F(t) — 
F p (t) — e - 2Rcr ( t ). if the single-particle r(t) may be dy- 
namically suppressed, i.e. if the spectrally shifted bath 
response Gjj(u>j + Aj) is small enough, this F(t) will be 
kept close to 1 (Fig. If IIP symmetry is imposed, the 
particles become effectively independent in terms of their 
coupling to the bath. Hence, collective coupling to the 
bath, which is a prerequisite for a DFS, is not formed in 
this case. In order to impose decoherence-free conditions, 
we extend the treatment to multilevel particles 0. We 
will show that the use of auxiliary levels and local mod- 
ulations can result in a decoherence-free singly- excited N 
qubit system, in a general and realistic coupling scenario. 
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FIG. 2: Fidelity as a function of time (in units of (lO^) -1 ): 
(a) overall fidelity, (b) correlation preservation, and (c) 
population preservation. (i) Global 7r-phase flips impose 
no symmetry (dashed) (Tj, n = l-lj^j.n/t = 1.0). (ii) inde- 
pendent identical particles (IIP) symmetry (dotted) (Tj, n = 
(0.75, 0.85, 0.95, 1.05), 3 >/7r = (0.834,0.806,0.836,0.82)); 
(iii) independent identical trapping (IIT) symme- 
try (solid) (tj,„ = (0.85, 0.85, 1.05, 1.05), = 
(0.924,0.9,0.945,0.91)). Top-right: Decoherence ma- 
trix elements at time t = 100. The parameters 
are: U\ = 0.5, ui2 = 0.6, 7 = 0.1, k r min = 1.0, 
t j>n = (0.7, 1.0, 1.06, 1.1), rj = (0.0, 0.1), 77„/7r = (0.0, 0.1). 

To this end we invoke three-level particles with ex- 
cited states |1) and \2). In a single three-level particle, 
an external field can impose an intraparticlc DFS if |1) 
and 1 2) decay at the same rate to the ground state |0); 
this DFS consists of |0) and the "dark" or trapping state 
|l("| , the anti-symmetric superposition of the two excited 
states l/v2(|l) — |2)). However, this intraparticle DFS 
would be destroyed in a system of N three-level particles 
that are coupled to the bath, and/or experience cross- 
decoherence. In order to remedy this problem, let us 
consider local modulations acting differently on the two 
excited levels within each particle. Such modulations can 
be tailored to impose what may be dubbed "independent 
identical trapping" (IIT) symmetry. This means that all 
particles acquire identical trapping states and become 
(effectively) independent, without cross-decoherence. As 
discussed above, N spectral shifts Aj determined by 
|etj(w)| 2 suffice to eliminate cross-decoherence, Eq. <|15ll . 
Under these conditions, N destructively interfering pairs 
of local fields will cause trapping in each particle (Fig.^J, 
nji). This would result in an A^-dimensional DFS, com- 
posed of N particles sharing single-excitation with anti- 
symmetrically superposed excited states. Under the IIT 
symmetry the decoherence matrix is block diagonal, each 
block corresponding to particle j with levels n, n': 

4w(*)=%'^(*) (17) 
Both disentanglement and population loss are nearly 
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completely eliminated within the multipartite subspace, 
resulting in F(t) very close to 1 (Fig. 

As an example, the IIT recipe will be analyzed for 
two three-level particles, where each level of each parti- 
cle experiences different coupling to the bath, and cross- 
decoherence exists. Let us take the initial state of par- 
tides A, B to be |* s (0)) - (\-) A \0) B ± \0) A \-) B ) /v% 
where (j = A, B) = (jl^ - \2)^ jyfr. The goal will 

be to keep this state intact, |\I>s(i)) ~ |^s(0)), by pre- 
venting cross-decoherence and imposing intraparticle de- 
structive interference. The bath response matrix will be 

taken to be $jj>, nn ' (*) = ld n d n > - — fe^f~+FT) ~ wnere 
7 is the single-particle unperturbed decay rate (satisfy- 
ing Fermi's Golden Rule [13): d n = cos?/„, with r) n being 
the angle of transition dipole, tj <n is the correlation time 
of level n of particle j, kg is the wavevector of the reso- 
nant transition |0) 1 1(2)), r min is the distance below 
which the particles are identically coupled to the bath 
and rjji = \rj — Tji\, where rj is the position of parti- 
cle j. This model may describe the distance- dependent 
cross-decoherence of either radiatively or vibrationally re- 
laxing atoms or ions at different sites in a trap, lattice or 
cavity [p3 , IT^| . The parameters of the pulse sequence in 
Eq. © for creating the IIP or IIT symmetries were cho- 
sen such that the faster decaying particle experienced the 
stronger and/or faster pulses, i.e. Oj,n/Tj,n were adapted 
to the bath response parameters ^jj, nn of each particle. 
The pulse sequences were chosen to be different enough 
for each particle in order to eliminate the desired cross- 
decoherence terms, as per Eq. (|15f) . Individual pulses 
were selected to obey the other IIP or IIT requirements 
described above. The very high fidelity achievable by IIT 
is seen in Fig. [21 



We may apply these considerations, for example, to 
cold ions or atoms in a cavity [13J , whose radiative deco- 
herence rates depend on their position within the cavity 
and may vary by as much as 15%. This implies that 
placing several ions in the same cavity breaks their mul- 
tipartite symmetry. AC Stark shifts are then an effective 
tool to restore the symmetry. To impose the IIP symme- 
try, one requires an impulsive phase modulation at a rate 
exceeding the cavity-mode linewidth, 1/r > T (typically 
~ YlGBz). In order to impose the IIT symmetry, one 
can either use two excited states of the ion, or two ions 
for each qubit, such that their singlet Bell state and the 
ground state form the DFS Q . The modulations should 
be local, addressing each level and each ion separately, 
and adjusting the rate and phase of the impulsive phase 
modulation so as to impose the IIT multipartite symme- 
try, thus creating an iV-qubit DFS. Recent experiments 
13] indicate that such addressability is feasible either by 
spatial or spectral resolution. 

To conclude, we have shown that local modulations 
are generally far more apt than identical (global) mod- 
ulation to impose multipartite symmetry on an other- 
wise completely asymmetrically relaxing system. In par- 
ticular, local modulations can impose the IIT symme- 
try, in a system of TV three-level particles, thus creat- 
ing a decoherence-free singly- excited N -qubit system. The 
general formalism presented here spans the entire range 
of possible coupling scenarios and modulation schemes. 
These results imply that the challenge of multipartite de- 
coherence, plaguing quantum information transmission 
and storage, may be successfully met in a variety of ex- 
perimental situations. 
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